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Abstract 

We present a self-contained analysis of theories of discrete 2D gravity coupled to 
matter, using geometric methods to derive equations for generating functions in terms 
of free (noncommuting) variables. For the class of discrete gravity theories which corre- 
spond to matrix models, our method is a generalization of the technique of Schwinger- 
Dyson equations and is closely related to recent work describing the master field in 
terms of noncommuting variables; the important differences are that we derive a single 
equation for the generating function using purely graphical arguments, and that the 
approach is applicable to a broader class of theories than those described by matrix 
models. Several example applications are given here, including theories of gravity cou- 
pled to a single Ising spin (c = 1/2), multiple Ising spins {c = k/2), & general class of 
two-matrix models which includes the Ising theory and its dual, the three-state Potts 
model, and a dually weighted graph model which does not admit a simple description 
in terms of matrix models. 
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1 Introduction 



A major goal of theoretical physics is the construction of a self-consistent quantum theory of 
gravity coupled to matter in four dimensions. Because such a theory has proven to be very 
difficult to construct, it is interesting to consider toy models of gravity coupled to matter in 
lower dimensions, in the hope that such simplified theories can provide information about 
the nature of quantum gravity. In particular, it has been possible to construct consistent 
quantum theories of (Euclidean) gravity in two dimensions coupled to conformal matter. 
These theories have been the subject of considerable interest over the last decade, partly 
because they are the simplest available theories of gravity coupled to matter and partly 
because of the relevance of these models to string theory. Two dimensional gravity theories 
of this type can be defined either through the continuum Liouville theory approach 0, ||, 
or through discrete models of dynamically triangulated gravity (for reviews see @, i|, 

The study of discretized two-dimensional gravity theories has a long history §. A 
major breakthrough occurred when it was realized that a summation over "triangulations" 
of a 2-manifold by equilateral polygons (not necessarily triangles) can be described by a 
zero-dimensional field theory, or matrix model [9-14]. By writing a discretized 2D gravity 
theory as an integral over hermitian matrices, it becomes possible to use analytic methods to 
calculate correlation functions in these models, and to show that in an appropriate continuum 
limit these models appear to correspond to continuum theories of Liouville gravity coupled 



to conformal matter field (for a review see ||T5|). The matrix model technology which was 



developed also made it possible to sum the partition function for these models over Riemann 
surfaces of all genera, giving nonperturbative results for many theories of this type [IIB, |18 



One approach which has been used to analyze matrix model theories involves the use 
of Schwinger-Dyson equations |20|, These equations, derived by demanding the 



invariance of the matrix model correlation functions under infinitesimal changes of vari- 
ables, may be interpreted geometrically as "loop equations" describing the deformation of 
two-dimensional geometries. The Schwinger-Dyson/loop equations are a powerful tool for 
analyzing the algebraic structure of correlation functions, although it has proven difficult to 
find solutions in more complex theories. Recently, the Schwinger-Dyson equations for matrix 
models have been used to discuss c < 1 string field theories [22-31]. 

In this paper we generalize the Schwinger-Dyson technique, developing a systematic 
method for deriving a single equation satisfied by the generating function in a very general 
class of discrete 2D gravity theories. Our approach is based on the combinatorial methods 



pioneered by Tutte [32], who found an algebraic description of one class of triangulations by 



writing a discrete recursion relation for the number of triangulations of a 2-manifold with 



a fixed number of boundary edges (see also [33, M). We show how this type of recursion 



relation can be derived in the presence of discrete data defined on the triangulation, rep- 
resenting matter fields in the continuum. Using this technique we are able to describe a 
variety of theories, including theories such as the recently described dually weighted graph 
models [^S], ^ which do not always admit a simple description in terms of matrix models. 
For matrix model theories, this formalism can be used to describe generating functions for 
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correlation functions with a general class of boundary conditions, which can lead to results 
which cannot be easily obtained through the matrix model formalism. 

In theories with matter, the generating function for the disk amplitude is a formal sum 
over all possible matter configurations on the boundary of the triangulations. We encode 
boundary data for the matter fields as a string of free (noncommuting) variables associated 
with the matter configuration on the edges or vertices of the boundary. In this way the 
generating function is an element of the free algebra generated by the variables representing 
different matter states. One of the principal results presented in this paper is the derivation 
of a "generating equation," satisfied by the generating function. This generating equation 
is closely related to the master equation which has recently been discussed by a variety of 
authors [37-42]. The generating equation describes the effect of removing a marked edge 
from a triangulation, which results in one of two basic geometric "moves" - the removal of 
the triangle attached to the edge or, if the edge is connected directly to another boundary 
edge, the removal of that pair of boundary edges. In fact, the correspondence between these 
moves and the generating equation is sufficiently precise that the equation may be written 
down directly upon consideration of the possible effects of the moves, without going through 
the intermediate step of an explicit recursion equation. 

Once we have a single generating equation in terms of noncommuting variables for the 
correlation functions in the theory, the next step is to try to find a subset of correlation 
functions for which a closed system of equations can be written in terms of commuting 



variables. This is analogous to the technique used in Ref. |^3|] to find the cubic equation 
satisfied by the generating function of homogeneous correlation functions in the c = 1/2 
matrix model which describes an Ising matter field coupled to 2D gravity; a description of this 



calculation in terms of noncommuting variables was given in Ref. |^2[. From a polynomial 
equation of this type it is straightforward if somewhat tedious to derive information about 
the scaling behaviour of correlation functions in the continuum limit. 

A principal focus of this paper is the development of a systematic approach to finding 
and solving closed sets of equations in commuting variables starting from the noncommuting 
generating equation. We give several examples of polynomial equations which can be derived 
in this fashion. In particular, we explicitly solve a 1-parameter family of two-matrix models 
which may be thought of as the Ising model coupled to gravity with a boundary magnetic 
field. This family of models includes both the Ising theory (in which Ising spins are located 
on the faces of the triangles) and its dual model (in which the spins are located on the 
vertices of the triangulation). The solution reproduces the cubic equations for the Ising and 
dual theories, and shows that the intermediate models satisfy a quartic. We also describe 
the generating equations for the 3-state Potts model and the 2'^-matrix model with c = k/2 
which describes k Ising spins coupled to 2D gravity; in both of these examples we are unable 
to find a closed system of equations in commuting variables, and we discuss the obstacles 
encountered in this approach. 

We do not discuss the continuum limit of any of the polynomial equations that are 
derived in this paper. All the equations are at most quartic in a single commuting generating 
function, and contain a finite number of functions that correspond to boundary data for the 
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generating function. The continuum limit is given by the expansion of the parameters of the 
theory around their critical values, as governed by the appropriate polynomial equation; for 



an example see p2[. The continuum limit of some of the more interesting models considered 



in this paper will be discussed in Refs. [44, 45 



The approach we have taken in this paper is to develop the general formalism in a 
piecewise fashion, introducing each additional complication in the context of a particular 
model where extra structure is present. The outline of the paper is as follows: in Section 
2 we describe the geometric approach used throughout the paper by examining in detail 
the generating equation for pure triangulated 2D gravity, which can be written in terms 
of a single variable that encodes boundary information. In Section 3, we consider gravity 
coupled to a single Ising spin, and associate boundary configurations with words in a free 
algebra generated by noncommuting variables. We derive the cubic equation governing 
homogeneous boundary configurations in this model. We also give a generating equation 
for a more general class of boundary variables in this model which incorporates information 
about nearest-neighbor spin pairs. In Section 4 we study the Ising theory on the dual lattice, 
and derive the cubic equation for the homogeneous disk amplitude in this formulation of the 
model. In Section 5 we find the generating equation for a general two-matrix model with 
cubic interactions, which includes the Ising theory, dual Ising theory, and the 1-parameter 
family of models describing the Ising theory with boundary magnetic field. We also outline 
the derivation of a quartic equation for the homogeneous disk amplitude in these models. 
Some of the results mentioned in this section will be described in greater detail in a separate 



paper |^5[. In Section 6 we examine the 3-state Potts model, which serves as an example of 
the obstacles faced by this kind of approach. In Section 7, we describe the generating equation 
for the theory of triangulated gravity coupled to multiple Ising spins {c = k/2). In Section 
8, we consider a more general class of discrete gravity models, and derive the generating 
equation for a simple dually weighted graph model. Section 9 contains concluding remarks. 



2 Pure Gravity 

In this section we consider pure 2D gravity in the absence of matter fields; this corresponds 
to the generic one-matrix model. We use this simple model to present an introduction to 
the geometrical techniques that are applied to solve more complex models in subsequent 
sections. We derive the well known equation for the disk amplitude generating function, in 
the case of cubic interactions fl^ . Although the generating function in this model is written 
in terms of a single commuting variable, we show how this result may be written in a form 
that generalizes to free variables. Finally, we derive a more elaborate system of generating 
equations for higher-genus amplitudes. 
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2.1 Review of model 



The partition function for pure 2D gravity is formally defined as 



r y DiffM / 

where the genus h of the orientable 2-manifold Ai is summed over and 

Se[9,u] = / d'x^{-R + 2 A) . (2) 

A discretized version of this theory can be constructed by summing over all triangulated 
2-manifolds. This discrete theory has a partition function 

where g and N are coupling constants that replace A and G, x(^) is the Euler character of 
the triangulation A, 5* (A) is the symmetry factor of A, and ?t,(A) is the number of triangles 
in A. The matrix model expression for this partition function is given by 



Z = JdU exp 



ixr f/2 - ^Tr f/^ 
2 3 



(4) 



where U is an hermitian N x N matrix. 

We are interested in computing the generating function ^{u,g) for the disk amplitude, 
defined by 

oo oo 

g) = Y: PkigV = E ^(k; n)g-u' , (5) 

fc=0 k,n=0 

where pk is equal to the partition function summed over all triangulations having k boundary 
edges. The correspondence with the matrix model language is that A/'(fc; n) is equal to the 
number of planar diagrams (up to symmetry factors) with n trivalent vertices and k external 
legs; these may be thought of as the dual graphs to the triangulations. Thus, pk is given by 
the matrix model expectation value 

Pk = (Tt {U')) . (6) 

so that 

(Note: in this paper all matrix model expectation values are evaluated in the large A^ limit.) 

There are two simple ways to solve for the function ^[u,g). In Ref. [^, a recursion 
relation is derived for the number of triangulations in which no two edges join the same pair 
of vertices, no internal edge is connected to two external vertices, and no two triangles have 
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three vertices in common. From the point of view of the dual graphs, this is equivalent to 
demanding that the graphs be built from connected, one-particle irreducible diagrams with 
no self-energy corrections to the propagator. Similar recursion relations are discussed in Ref. 



3^. Alternatively, ^{u,g) can be evaluated by computing the functional integral over 
matrices. In this case, all Feynman graphs are summed over automatically, and the coeffi- 
cients of ^{u,g) can be obtained by saddle point methods or by the method of orthogonal 
polynomials |10|- Although the numerical values of the coefficients Af{k;n) are different 
in the Tutte and matrix model schemes, the universal behaviour of ^{u,g) (which governs 
the continuum limit) is the same. We shall follow Tutte in deriving recursion equations for 
the number of triangulations, but we adopt the matrix model conventions for counting: our 
dual graphs need not be connected nor one-particle irreducible, and we allow self-energy 
corrections to the propagator. 



2.2 Generating equation 

We now proceed to derive a recursion relation for the coefficients //{k; n), which leads to a 
single generating equation satisfied by g). We consider a triangulation of a disk D with 
n triangles and k external edges, one of which is marked. The given triangulation can be 
disassembled step by step, by repeatedly removing the marked edge and marking a new edge 
on the resulting, smaller triangulation. (A similar approach has been used recently to derive 
a discrete Hamiltonian operator for string field theory in the proper time gauge |2^, ^ |47||.) 
Each time a marked edge is removed, this has the effect of either removing a triangle, or 
removing a pair of boundary edges. In this way, for each triangulation there is a unique 
sequence of "moves" reducing to the trivial triangulation, a single dot. By considering all 
possible sequences of moves, we can recursively count the number of triangulations J\f{k; n). 

Some conventions are necessary to make the recursion relations well-defined. When the 
marked edge is a boundary of a triangle, one removes the triangle attached to the edge; the 
resulting triangulation has one less triangle, but one additional edge. We use a convention 
that marks the edge which was counterclockwise from the original marked edge in the triangle 
that was removed (see Fig. [^). In the case where the marked boundary edge is connected 
to another boundary edge, one should remove the pair of edges, leaving two triangulations. 
On each of the remaining triangulations, we choose the convention of marking the edge next 
to edge which was just removed (see Fig. |^). In Fig. ^ we illustrate the results that can be 
obtained upon removing various possible marked edges from a simple triangulation. 

The number of triangulations with n triangles and k boundary edges is therefore equal 
to the number with one fewer triangle and one additional edge, plus the sum of all possible 
combinations of two triangulations with the same number of triangles but a total of two 
fewer boundary edges: 

k-2 n 

M{k] n) = M{k + 1; n - 1) + ^ ^ U{i; n')U{k -2-i-n-n') . (8) 

i=0 n'=0 

In order that this recursion relation completely determine M{k; n), it is necessary to impose 
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Figure 1: The marked edge is the boundary of a triangle. The triangle 
is removed. Of the two new edges, the one further counterclockwise is 
marked in the new triangulation. 




Figure 2: When the marked edge is connected to another boundary edge, 
the triangulation is split into two disconnected pieces. Each new tri- 
angulation has a marked edge next to where the removed edges were 
attached. 



boundary conditions: there is one trivial triangulation, so that A/'(0, 0) = 1, and in addition 
we specify that Af{0;n) = for all n > 1. These conditions suffice to recover the counting 
of the matrix model. 

Using (^), it is straightforward to rewrite (|^) as an equation for the generating function 

g) = l + g) - pi{g)u - 1) + u^<l>\u, g) (9) 

u 

The first term on the right hand side of @ fixes the boundary condition A/'(0; n) = 6n,o- The 
second term corresponds to the first term in (|^), and represents the removal of a triangle, 
as in Fig. |l|; the negative terms are subtracted from ^{u,g) since no move that removes a 
triangle can result in a triangulation with fewer than two edges. The third term corresponds 
to the second term in (H). 

It is important to note that (P) cannot be solved for $ directly, since it contains the 
function pi{g) which is a boundary condition for ^{u,g) that must be computed separately. 
Equations of this general form have been studied extensively in the literature, and it is 
possible to show that the requirement that g) have a well defined power series expansion 
in u and g around the origin is sufficient to completely determine the power series pi{g) from 
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Figure 3: A simple triangulation is shown at top, with numbers repre- 
senting different edges which could be marked. The six diagrams below 
are the results of removing the appropriate marked edge, and show which 
edge of the new triangulation becomes marked. 
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(y). We refer the reader to Refs. ^ for details of this procedure for the present 
example, and for the more complicated cases that we shall encounter in later sections. 



2.3 Graphical interpretation 

The generating equation (^) may be written more compactly by introducing a special deriva- 
tive operator, which will turn out to be useful in more complicated matrix models where $ 
will be a function of free variables. We define the operation of du on a polynomial in u so 
that it removes one power of u from terms proportional to (without introducing a factor 
of k), and annihilates terms which are independent of u Thus, 

oo oo 

duY.Cku'' = Y.CkU^~^ . (10) 

k=0 k=l 

Such an operator does not obey the usual Leibniz rule; rather, 

9„($^) = (9„<l>)^ + ^{u = 0)9„^ . (11) 
Using this definition, we have 

5„$ = -($-l), 9^$ = 1($-1-Mpi) , (12) 
u 

and d^) reduces to 

$(«, g) = l+ gudl^iu, g) + u^e{u, g) . (13) 

In this notation, the correspondence between the terms in the generating equation and 
the graphical representation of the basic moves becomes precise, in the sense that every 



symbol appearing in eq. (|T3D has a direct interpretation. On the right hand side there are 
three terms, corresponding to a boundary condition, the move which removes a triangle, 
and the move which removes two identified boundary edges; each of these three types of 
terms will appear in the more elaborate generating equations we present in the remainder of 
the paper. The right hand side of each such equation will begin with a one, which sets the 
number of triangulations with no boundary edges equal to one. The terms corresponding to 
removal of a triangle will have a factor of the weight for a single triangle (in this case, g) and 
an operation in free variables (in this case, udf^ which represents replacing the single marked 
edge with two new edges.Q The terms corresponding to removing two identified boundary 
edges will contain a variable for each boundary edge (in this case, two factors of u) and two 
functions (in this case, two factors of $) representing the two triangulations that remain 
after the removal of the boundary edges. In more complicated models these terms may be 
multiplied by coupling constants representing the interaction between matter fields, as we 
shall see in the next section. 



^The notation ud^ may seem to suggest replacing two edges with a single one, rather than the other way 
around. This is merely because we have been using language in which the triangulations represented by $ 
are disassembled step by step, rather than being built up from nothing; but either sense is equally legitimate. 
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Eq. ( [TB| ) is closely related to the Schwinger-Dyson equation for the 1-matrix model, which 
can be written as 

du^ = + gdl^ . (14) 



It is easy to see that acting by du on (|T3]) yields (|T^). The generating equation includes the 
information contained within the Schwinger-Dyson equation, as well as additional boundary 
conditions. This will continue to be the case for more complicated theories described by 
multi-matrix models, where we shall derive a single equation equivalent to the entire set of 
Schwinger-Dyson equations plus boundary conditions. 



2.4 Higher genus 

It is reasonably straightforward to extend the recursion relation (^ from the disk amplitude 
to an amplitude for a higher genus surface with an arbitrary number of boundaries. We 
imagine that the set of boundaries is ordered, and that each one has a marked edge, and 
derive a recursion equation by removing the marked edge on the first boundary loop. The 
number of triangulations of a surface with genus G and with b boundaries, represented by 

Af^'\h,k2,...,h;n) , (15) 

with all ka > 0, is given by the equation 

Xg\k,,k2,...,k,;n) = J\f^'\k, + l,k2,...,h;n-l) 

b-l G fci-2 n 

0=0 {a-i} H=0 i=0 771=0 
b 

+ Y kaJ^G'^\ki + ka-2,k2,..., fca-1, ^a+l, ■ ■ ■ , h] U - 1) 
a=2 

+ E^G-i'^(^'^i-^-2,A;2,...,A;„;n) 

i=0 

(16) 

where {(jj} in the second term represents the set of all permutations of (2, 3, ... , b) such that 
0"! < < . . . < (Tfl and (Ja+i < o"a+2 < • • • < and we set NQ\ki, k2, . . . , kf,; n) = 

when any of the kaS are zero (except that, as before, J\fQ^\0; 0) = 0). 

Once again, each term in this expression has a geometric interpretation (see Fig^. The 
first term on the right hand side corresponds to the removal of a single triangle (Fig ^a), 
while the remaining three correspond to the removal of an identified pair of boundary edges. 
The second term is related to the second term in (|]), derived from the splitting of a single 
triangulation into all possible combinations of two (Fig Qd). The third term corresponds 
to the removal of two identified edges which belong to different boundary components; the 
result is to decrease the number of boundaries by one (Fig |^c). The last term comes from 
the removal of two identified edges belonging to the same boundary component, such that 
the resulting triangulation does not split into two; the result is to increase the number of 
boundaries by one, while decreasing the genus by one (Fig ^). 
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(a) 



(b) 




(c) 




Figure 4: For a higher genus surface, there are two new moves (c) and 
(d), as well as the two moves (a) and (h) needed for the disk amplitude. 



One can define generating functions 

oo oo 
G=0 {kj},n=Q 

for triangulations of a surface with h boundaries of arbitrary genus. It is now possible to 
convert the recurrence relation (0) into an equation that relates amplitudes for surfaces of 
arbitrary genus with different numbers of boundaries: 

^^^\ui,...,Ub,g,h) = 

6bi+guidl^<l>^^\ui,...,Ub,g,h) 

b-l 

E E ^^""^^^ui, u^„..., u^^,g, h)^^''~''\ui, u^^^„ ...,Ub,g, h) 
^ d V 

+ E"il~ uiUjDu^^uj^^^''^\ui, . . . ,Uj-i,Uj+i, . . . ,Ub,g,h) (18) 

b 

+hul E "^^^^^^ (mi, • • • , Ml, Uj, ...,Ub,g,h) , 
i=2 

where {cTj} is defined as before, the operator D^^ ^ is defined by 



^) (mi, . . . , Mj-i, Mj+i, ...,Ub,g, h) 



^) (mi, . . . , ...,Ub,g,h)- 1) (m2, • • • , ^6, 5', ^) 



Ui — Uj 



(19) 
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and d/dui is an ordinary derivative operator. 

Equation (|T8|) is useful in part because it contains a set of linear relations that permit an 
iterative solution for generating functions of arbitrary genus and with an arbitrary number 
of surfaces starting from the solution of Eq. (^) for the disk. For example, the simplest 
relation contained in that equation is obtained by projecting out the genus zero component 
and taking 6 = 2. This gives a relation between the disk and cylinder generating functions 
(the components $o^^ and $o^^ of and $(^)): 

d (ulu2^\^\u,,g)-ulu,^'i\u2,g)\ (20) 



U2 



dU2 \ Ui — U2 



Another simple relation exists between the generating functions for the disk, the cylinder 
and the single-boundary torus (the component taking 6 = 1): 

<l>i'HMi,(7) = 2M2$«(Mi,(7)$S'Hwi,^?)+^?Wi5ni'^'i'H^ • (21) 

This procedure may be continued to obtain equations for arbitrary genus and number of 
boundary components. 

3 Ising matter 

We now consider 2D gravity coupled to the Ising model. This model corresponds to the 
c = 1/2 minimal model CFT coupled to Liouville gravity, and has a simple description as a 
two-matrix model. As in the previous section, most of the results we describe here have long 
been understood in terms of the matrix model description of the theory; the novelty of our 
presentation lies in the purely geometric derivation of a single generating equation written 
in terms of noncommuting variables. 

3.1 Review of model 

The coupling of the Ising model to gravity is described by summing over all triangulations 
of a closed 2-manifold just as in the pure gravity model. Now, however, each triangle carries 
a label U or V , corresponding to Ising spins up or down respectively. The weight of a given 
triangulation A is defined to be 

1 / 1 \ mm(A)+ot(A) / ^ \ Mu(A) 

wiA) - ^A'-V- ijhr) (l^) 

where the notation is the same as in (^, with the addition of a new coupling constant c; 
M«(A) is the number of edges in A separating two U triangles, and similarly for uv{/S) and 
(A). The couphng c is related to the usual Ising model coupling J by 

c = . (23) 
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We use c rather than J because this is the natural notation in the description of the theory 
as a 2-matrix model, with partition function given by |]48| , |49 | 



Z = Y, = [dUDV exp {~NS{U, V)) 



(24) 



with 



S{U,V) 

where U and V are N x N hermitian matrices 



-Tr f/2 + -Tt V'^ - cTr UV - -(Tr + Tr V^] 
2 2 3^ ^ 



(25) 



3.2 Generating equation 

We would like to analyze the generating function associated with a disk, as in the pure 
gravity theory discussed in Section 2. However, in this model it is necessary to specify the 
configuration of Ising spins on the boundary rather than just the number of edges. Thus, 
we must associate with each boundary a length /, and a string of / labels, each of which is 
either a u or a f . These labels may be thought of as spins lying outside the boundary of the 
triangulation, meaning that each boundary edge contributes a weight 1/(1 — c^) or c/ (1 — c^) 
depending on whether the spin on the boundary edge is the same as or different from the 
spin on the triangle adjacent to that edge; although this may seem to be a counterintuitive 
convention, it will simplify the following analysis. By choosing an orientation and a marked 
edge on the boundary, we have a natural ordering to the boundary labels; this ordering is 
given by beginning with the marked edge and proceeding around the boundary. It is now 
natural to consider u and v as noncommuting variables; the generating function for the disk 
can then be written as (throughout this section we suppress the dependence of $ on the 
couplings g, c and the noncommuting variables u, v) 

$ = 1 + Puiu + V) + Puuiu"^ + v"^) + PuviuV + VU) 

(26) 

+Puuu{u^ + V^) + Puuv{uUV + UVU + VUU + UVV + VUV + VVu) + . . . 

The coefficients Pw{u,v) represent the disk amplitude for triangulations with boundary con- 
dition specified by w(m, f ), a word in the free algebra generated by u and v {i.e., an ordered 
product of m's and f 's). Notice that, due to the symmetry between plus and minus spins, 
the expectation values are related by Pw{u,v) = Pw{v,u)- (In the models we shall study in 
subsequent sections, this will no longer be the case.) There are also symmetries under cyclic 
permutation and reversal of w{u,v). In matrix model language, Pu}{u,v) is just the usual 
expectation value of a trace of a product of hermitian matrices (in the large limit): 

Pn,{u,v) = 4(Tr^(f/,^)) , (27) 
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Figure 5: A diagram contributing a factor of Ng^c^/{1 — c^)^ to the 
expectation value (Tr U^) = Npuuuu- 



and $ is the generating function of planar Green's functions WA HI- 

1 oo 

* = E (Tr («f/ + 

1 (28) 
= nT. w{u,v)(Tiw{U,V)) . 

w{u,v) 

For example, Figure ^illustrates a diagram contributing a factor of Ng'^c^ /{I — c^)^ to the 
expectation value (Tr If^) = Npuuuu- 

We can now, as in the previous section, use purely geometrical arguments to derive a single 
generating equation satisfied by taking advantage of the interpretation discussed after 
Eq. ([13| ) to write down the generating equation directly, without first deriving a recursion 
relation. We have associated the first symbol in a word w{u, v) with a marked edge on 
the boundary of a triangulation. If we remove this edge, then as in the pure gravity case 
there are two possibilities: either the edge is identified with another edge of the boundary 
(which must also be removed), or the edge is attached to a triangle which can be removed. 
Unlike the pure gravity case, however, in each of these cases there are multiple possible spin 
configurations which must be included. 

It is useful to define derivative operators du and that are similar to the derivative 
operator defined in the previous section. (Operators of this form were introduced in Refs. 
1^, 0.) If du (dy) acts on a word with a u (v) at the leftmost position, that u (v) is removed; 
if a w {u) is at the leftmost then the term is annihilated. For example: 

du{uuuvv + vvu) = uuvv, dy{uuuvv + vvu) = vu . (29) 

In terms of these operators, the generating equation is 



$ = 1 + — 



+ f $f $ + cM$f $ + cf $u$ + gudl^ + gvdl^ + cgvdl^ + cgud^^ . 

(30) 
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Figure 6: The removal of a pair of edges of opposite spin results in a 
term cf on the right hand side of the generating equation. The 
removal of a triangle with a spin at its center not equal to the spin on 
its edge results in a term cgudl^. Note that the two new edges carry 
the spin of the triangle that was removed. 



Again, each term corresponds to a possible outcome of removing a boundary edge. The first 
four terms in the square brackets are associated with the removal of two identified edges: a 
marked u identified with another m, a marked v identified with another w , etc. (Note the 
factor of c appearing when a m is identified with a. v.) The other four terms correspond to 
removing a triangle (and therefore have a factor of g)\ they represent the four possibilities 
of a M or f marked edge being identified with a m or f triangle. A diagrammatic description 
of the fourth and eighth terms inside the brackets is given in Figure ^. 

The generating equation (|30|) is closely related to the Schwinger-Dyson equations for the 



corresponding matrix model. The Schwinger-Dyson equations are given by 

a„$ = $M<I> + gdl^ + cd^^ 

(31) 



It is fairly easy to check that these equations are equivalent to ( POD except that they do not 
determine the leading constant in $. Thus, the single generating equation contains all the 
information contained in the pair of Schwinger-Dyson equations. Having a single equation for 
$ in terms of noncommuting variables simplifies the subsequent analysis of such a system, 
as we shall see in the following sections. Finally, note that it is straightforward to write 
down the generalization of the higher genus generating equation (p!8D in the presence of Ising 
matter. 
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3.3 Solving the generating equation 

Now that we have a single generating equation describing the set of all disk amplitudes, we 
would like to find a closed form solution for $ in terms of noncommuting variables. Unfortu- 
nately, however, finding a complete solution to such a nonlinear equation in noncommuting 
variables seems to be a difficult problem It is possible to directly compute any 

specific disk amplitude with fixed boundary length / and a finite power of g by repeated 
application of (|30|). However, to understand the continuum limit of the theory we need to 
have information about amplitudes with arbitrarily large boundaries. In certain cases it has 
been found possible to find a closed set of loop equations for matrix models which allow the 
computation of generating functions with certain constraints on boundary configurations 
|]50| , pT| , We shall now demonstrate in the generating function approach a simple sys- 
tematic procedure by which such closed systems of equations can be found and solved. The 
following analysis closely parallels that of |^ . 



By factoring out leading terms which depend on f , the generating function $ can be 
expanded in terms of functions (pwiu^v) which only depend on u: 

v) = (j){u) + v(j)^{u) + v'^(j).,y{u) + uvdu4>v{u) 

{o2) 

+v^4>vvv{u) + vuv(f)vuv{u) + uvvdu(f)vv{u) + uuvdl(j)^{u) H 

The functions (t)w{u,v){u) are defined by 

<Pwiu,v) = [widu,dy)^]^^Q , (33) 

where w{u, v) begins and ends with v. (If w{u, v) begins with a u, we can use the cyclic sym- 
metry of the generating function to replace the corresponding function 4>w(u,v) by a derivative 
of a lower-order function, as we have done in (|32D .) Thus, w{u,v)(f)w{u,v) is that piece of $ 
containing words consisting of w{u,v) multiplied on the right by a string of m's. (Note that 
(0) actually extracts the term with the m's and f 's in the opposite order, but the functions 
are equal by symmetry.) 

Considering only those terms on each side of ( pOD with a fixed initial string followed by 
m's, we can write a closed system of equations for a subset of the functions (pw{u,v){u). For 
example, taking only terms with u's, we have 

= 1 + YZ~^ ("V^ + 9u9l(f) + cgu(f)^^) . (34) 
Taking terms of the form vu^, uvu\ v'^u^ give respectively 

du(l)v = ^_^2 ("00^' + c0 + gdl(l)v + cg(f)^^^^ (35) 

(pvv = z ^ (0 + CM0^0 + g(j)y^^ + cgdlq 

1 — ^ 
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This gives us 4 equations in 4 unknowns, each of which is a function of the single commuting 
variable u. These 4 equations can be combined to give a single equation for and its 
derivatives: 

= gu^(f)^ + (cM^ — gu + g^)4P' — g"^ + {c^u — cu — gu'^(j))du4> 

(36) 

+ {gu + cgu + g^u^(j))dl(f) - 2g^udl(j) + g^udf^(f) . 

As was done in Eq. (|T^) for pure gravity, we can write the derivatives of (j) in terms of (j), 
u, and a finite number of the individual coefficients Puu, and so on. This allows us to 
rewrite (|36| ) as a pure cubic in 0: 

/30' + /20' + /l0 + /o = O , (37) 

where 

/o = -g^ + g'^{2-gpi)u-g{l + c-2gpi + g^P2)u'^ 

- {g^ + g{l + c)pi - 2g^p2 + g^Ps - c(l - c^)) 
h = g-^-2g'u + g{l + c)u'-{c{l-c')+g')u' + g{l-gp,)u' (38) 
/a = u^{2g^ -2gu + cu'^) 

h = gu' 

Notice that the coefficients Pu and Puu appearing here are not independent; the order u terms 



in (|34|) give 

Puu = — — —Pu ■ (39) 
9 



The cubic (^TD describing homogeneous (fixed spin) disk amplitudes for the Ising model 
coupled to gravity is well known and has been derived using a variety of matrix model 
methods |5^, 51, i3, 42 1. The derivation we have given is completely independent of the 
matrix model formalism, and depends only on the underlying geometric theory. The existence 
of a single equation for the noncommuting generating function ([30| ) rather than a pair of 
Schwinger-Dyson equations has simplified our analysis. 

It is straightforward given the disk amplitude to obtain expressions for generating func- 
tions in commuting variables for boundary data with multiple domains of like spins. One 
proceeds iteratively, deriving an expression for the two domain amplitude in terms of the 



homogeneous amplitudes and so on |^3[. Defining t(u,v) as the two-domain generating 



function, including all terms homogeneous in u and v, it follows from ( pOD that 

T — (j)(v) = U(f>(u)uT(u,v) + CUT(u,v)v(j)(v) 

1 — L 

g ^g^ 1 
H — {t{u, v) - U(t)^{v) - 0(f)) H 5-(r(M, v) - v(j)y{u) - 0(m)) (40) 

which is linear in t{u,v). (Note that in this equation u and v act as ordinary commuting 
variables, and we have converted the derivative operators into their algebraic expressions.) 
Similar relations can easily be derived for the higher domain generating functions. 
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3.4 Alternate boundary conditions 

We have shown how noncommuting variables can be used to encode the boundary data used 
in a sum over triangulations coupled to Ising matter. However, the decision to describe the 
boundary as a sequence of plus and minus spins is not unique. In this subsection we briefly 
describe an alternative description of the Ising theory, in which the boundary condition 
is specified by labelling the vertices between boundary edges as either "stick" or "flip," 
depending on whether the spins on either side are like or unlike. We outline the derivation 
of the cubic equation corresponding to all-stick boundary conditions, which is equivalent to 
(^). Although we do not derive any results here which cannot be found using standard 
matrix model methods, there are certain correlation functions which can be derived using 
this alternative formulation which appear to be inaccessible using the matrix model language. 
An example of such a calculation appears in p4| . 

The generating equation in the flip/stick variables may be derived using the kind of 
geometric argument used to find (^). For example, using the symbols / and s for flip and 
stick, the term cgf'^dg^/{l — c'^) will arise when a marked edge bounded on each side by a flip 
is connected to a triangle, which when removed reveals two edges (and thus three vertices) 
labelled by sss. This corresponds in the u, v variables to removing a plus spin which is 
surrounded by two minuses and connected to a minus triangle, or the same operation with 
plus and minus interchanged. Taking all possible such moves into account, we obtain 



$ = 1 + 



1 -C2 



g[{s + cpds) + {cfs + sf)df]dl^ + g[{fs + csf)ds + {cs^ + P)df]dsdf^ 



+ 2 
1 

+ 2 



[cs^ + f){ds^f + df^s) + {fs + csf){2 + ds^S + df^f)\ (M$ + Sdf^) 
[s^ + cf){2 + ds^s + df^f) + {cfs + + df^s)] (2 + sds^ + fdf^) 



+g{l + c)sps 



(41) 



The leading term of $ is 1, but all other terms appear twice corresponding to the symmetry 
induced by interchanging pluses and minuses on the boundary. 

From here, we can define as before a set of functions (t>w{sj){s) and look for a system of 
equations which close. One such set of four equations is given by 



1 + - 
1 



-s^ (2 + sds(t)~f + g [s'^d^(j)^ + spss) + cgs {s(j)fsf + Ps 



l-c2 
1 



-c{2 + sds(t)-f + g(j)fsf + cgd^s^ 



-s (2 + sds 



+ Ps 



(42) 



+c (2 + sds(j)-) + eg {s(f)fssf + Psss) + 9 {sd'^s^'ff + Pfsf 



1-C2 



ds4>- + -c (2 + sds4>-) {ps + s(j)ff) 
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+9<t>fssf + cgdl(pff 



where the homogeneous disk amphtude with all sticks is represented by These equa- 

tions are exactly equivalent to the four equations in (|34D and (^) as can be seen by noting 
the relations 

0- = 20-l, S(j)ff + Ps = 2(f)y, S(j)fsf + Pss = '2(j)vv, S(f)fssf + Psss = '2(j)vvv ^^^^ 
Ps Pgg 2pj,^, Psss '^Pvvv 



where the factors of 2 arise because of the interchange symmetry. The cubic equation ob- 
tained by solving (^) for is therefore exactly equivalent to the equation for the ho- 
mogeneous disk amplitude 4>{u). 

Again, it is straightforward to extend the generating equation (^T]) to arbitrary genus. 
This extension turns out to be useful for computing the correlation functions of operators 



on the sphere 44 



4 The dual Ising theory 

One of the most remarkable properties of the Ising model on a fixed lattice is the Kramers- 
Wannier duality symmetry, which states that the partition function of the Ising model on a 
lattice L at a fixed temperature T is equal to the partition function on the dual lattice L at 
a temperature T' ~ 1/T. By putting spins on the vertices of a random triangulation, rather 
than on the triangles, it is possible to construct a dual version of the random- surface Ising 
theory discussed in the previous section. The resulting model can be written as a two-matrix 



model, and is a special case of a more general 0{n) loop model ||5^, |5^, |J, ^ ^Tj. In this 
section we consider this dual model, and show that the homogeneous disk amplitude of the 
dual model satisfies a cubic equation very similar to (|36|); an equivalent cubic was derived in 
||56| using different methods. This homogeneous disk amplitude for the dual model is precisely 
equal to the disk amplitude of the original Ising theory with free boundary conditions. 



4.1 Dual model 

In the dual model, we associate a spin with each polygon in the dual triangulation. The 
model can be described in terms of combinatorial data on the original triangulation by taking 
advantage of the fact that any dual spin configuration is determined (up to an overall sign) 
by a loop graph on the original triangulation - that is, by a set of edges which divide regions 
of opposite spin (see Figure 0). Each edge of the triangulation (or the corresponding edge 
of the dual triangulation) is labeled with a F if it separates different spins {i.e., if a domain 
boundary passes through it), and with an X if it connects identical spins. We then sum over 
all loop graphs on triangulations by summing over all triangulations and labelings such that 
the boundary edges of each triangle either have 3 X labels, or 1 X label and 2 Y labels. The 



19 



Figure 7: A sample triangulation contributing to the dual model. The 
dual graph, consisting of polygons which meet at trivalent vertices, is 
shown with doubled lines. Putting the spins on the dual graph rather 
than the original triangulation results in a different coupling between 
the matter and gravity in the two models. The edges labeled Y ( denoted 
by a thick line inside the double lines) define the boundaries between 
domains of like spin, and thus determine the matter configuration up 
to a sign. Note that a boundary condition with all X 's in the dual 
theory is equivalent to a homogeneous boundary configuration for the 
dual spins. 
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weight given to each triangulation is 



W{A) 



S{A) 



jYX(A)-n(A) 



x(A) 



1 + C 



?/(A) 



(44) 



where x{A), y{A) are the number of X,Y edges in A respectively. Putting the spins on 
the vertices of the triangulation leads to a different coupling of the Ising model to gravity 
than that considered in the previous section (see Figure We shall continue to refer to the 
model previously introduced as the Ising theory, and we will refer to the theory defined by 
Eq. ( ^4|) as the dual theory. 

A matrix model description of the partition function of the dual theory is |53, |56l 



Z = JdXDY exp (^-N 



C n (1 + C) o 

-Tr X^ + ^ ' Ti Y^ 



;(TrX^ + 3TrXF 




(45) 



It is a remarkable fact that the matrix model ( ^5]) is precisely equivalent to the original 
model (p^) [p6t This equivalence can be explicitly seen by replacing 

1 



X 



Y 



9 



(46) 



Because of this correspondence, the partition functions on a closed manifold for the two 
theories are precisely the same. Nevertheless, this does not by itself imply that Kramers- 
Wannier duality is preserved in the presence of gravity; the transformation (^) exchanges 
the disorder and spin operators, and takes homogeneous X boundary conditions in the dual 
model to free boundary conditions in the original model (and vice versa). A true duality 
would not exchange operators and would relate amplitudes with the same types of boundary 



conditions. In |44] we examine the issue of duality between these models more closely. 



4.2 Generating equation 

We again study the generating function for a disk. Boundary conditions on the disk are 
specified by a string of labels x or y, whereas in the Ising theory the labels are taken to live 
outside the boundary of the disk. Taking into account that there must be an even number 
of y^s on the boundary (since each y on the boundary is the endpoint of a loop), we can 
expand the generating function as 

$(x, y) = 1+ p^x + p^^x^ + pyyy"^ + Pxx^x^ + Pxyy{xyy + yxy + yyx) H , (47) 

where the coefficients may be thought of as expectation values: 

PK.,.) = 4(Tr^(X,r)) . (48) 
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We can now derive the generating equation for $ using the geometric techniques of the 
previous sections. We obtain 



$ = 1 



1 



1 + c 



gydyd^^ + gyd^dy^ + y^y^ 



(49) 



The terms in (^91) correspond to moves that are analogous to those for the Ising theory. Note, 
however, that the equation is not symmetric in x and y. Again there is a direct relation 
between ( ^91) and the Schwinger-Dyson equations 

^ \gdl^ + gdl^ + $x$ 



9.$ 



1 -c 

1 



gdyd^^ + gd^dy<^ + ^y^ 



(50) 



for the matrix model. 

As we did in the Ising theory, we may define functions which depend only on the single 
variable x: 

(p^(^,y){x) = [w{d^,dy)^y^^ , (51) 
where w{x,y) begins and ends with y. It is then possible to expand $(x,?/) as 



^{x, y) = (f){x) + yy(t)yy{x) + yxycpy^yix) 

+xyyd^4>yy{x) + yxxy(f)y^^y{x) + yyyy(j)yyyy{x) + 



(52) 



Bearing in mind the asymmetry between X and Y, we can use (^) to compute the 
generating function 0(x) for homogeneous X boundary conditions. As mentioned above, 
has two interpretations: it is identically equal to the generating function for free boundary 
conditions for the Ising theory, while it can also be interpreted as a generating function for 
homogeneous boundary conditions for the dual Ising theory of 

It is straightfoward to derive the closed set of five equations 



1 + 



1-c 



^yy 



^yxy 



^yxy 



d. 



xYyy 



1+C 
1 

TT'c 
1 

l-c 

1 

l-c 



(f) X + gxd^cj) + gx(j)yy 

(j) + 94>yxy + gdx(f)yy 
P^(j) + g(j)yxxy + gdx(j)yxy 
+ 
+ 



(53) 



^yxxy 



^yyyy 



x< 



'yy 



^yyyy 



+ 9d:. 



'x^yy 



from (^91). These equations can be combined to solve for 0. Once again replacing the 
derivatives of $ with algebraic expressions, we obtain a cubic equation for (p: 



/30' + /20' + /l0 + /o = O 



(54) 
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with coefficients 



/o = -4c^2(l + xp^ + x^p^^) + 2cgx{3 - c)(l + p^x) - 2cx'^{l - 

fi = Acg'^ - 2cgx{3 - c) + 2cx^(l - c^) + gx^{l - 3c) - 2g'^p^x^ 

f2 = x"^ [g"^ — gx{l — 5c) — 2cx'^{l + c)) 

h = gx^ 



g'^x^ 



(55) 



As in the case of the Ising theory, it is straightforward, given an expression for the ho- 
mogenenous disk amphtude, to derive expressions for amphtudes with nontrivial domain 
structure on the boundary through an iterative procedure. In this case, however, the objects 
of interest are the generating functions with a finite (even) number of F' s on the boundary. 
For example, the two-domain generating function t{xi^X2) is defined to be the generating 
function for all boundary configurations with a Y (the starting point of the domain wall) 
followed by a series of X's followed by another Y (the end point of the domain wall) followed 
by another series of X's. It follows from (EDI) that r(xi,X2) is given by the equation 



r(xi,X2) 



1 



1 + c 



{Xi)(j){x2)+gX^^{T{xi,X2)-(l)yy{x2))+gX2^{T{Xi,X2)-(l)yy{Xi)) . (56) 



An analysis of the continuum limit of the homogeneous and higher domain amplitudes for 



the dual model is given in Ref. [44| 



5 General cubic two-matrix model 

The previous two sections have described the Ising and dual Ising theories, which correspond 
in the continuum limit to two distinct descriptions of the c = 1/2 CFT coupled to quantum 
gravity. In this section we investigate the most general two-matrix model with cubic inter- 
actions. This general class of matrix models cannot be solved by saddle point or orthogonal 
polynomial methods, since the form of the interaction terms prevents the angular degrees 
of freedom from being integrated out using the methods of ES, However, by using the 



techniques of this paper it is possible to find the homogeneous disk amplitude for any model 
in this class. 

After a description of the general model and a discussion of the quartic equation satisfied 
by the homogeneous disk amplitude in the general case, we look at a simple one-parameter 
family of models which interpolate between the Ising and dual Ising theories. Every model 
in this family has a homogeneous disk amplitude which satisfies a quartic equation. These 
homogeneous disk amplitudes can be interpreted as the disk amplitude for the regular Ising 
theory coupled to gravity in the presence of a boundary magnetic field. This set of models is 
of considerable intrinsic interest, and we shall study it in more detail in a subsequent paper 
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5.1 Definition and generating equation 

The general two-matrix model with cubic interactions is specified by its action, which we 
may write in the form 

S = 7T7-r ^ Tr (aQ^ + 13R^ - 2-fQR) - J Ti^aQ^ + 3bQ^R + SdQR" + eR") . (57) 

The coefficients of the quadratic terms have been chosen to yield a simple form for the 
propagator, which in the {Q, R) basis is 



(58) 



■/3 7' 
,7 «. 

Note that a, which multiplies in the action, is what propagates R to R, and vice-versa. 

From the point of view of triangulations, this theory may be thought of as one in which 
the edges of each triangle are labeled either Q or R; weights are assigned to each triangle 
according to the labels on its three edges and to each pair of neighboring edges of triangles 
or boundaries. Thus, the weight of a given triangulation is 

W(^^'^ = ^_jVX(A)^ggg(A)^(jgr(A)^grr(A)grrr(A)^rr{A)^(jg(A)^gr(A) ^ ^gg^ 
S{A) 

where qqq{A) is the number of triangles with three Q edges, qq{A) is the number of times 
an Q edge is identified with another Q edge, and so on. The partition function is given by 

Z = Y, W^(A) = [bQBR exp [-NS{Q, R)] . (60) 

A 

As in the previous sections, we can use geometric arguments of the kind leading to ( ^0|) 



and (p9D in order to derive an equation satisified by the generating function associated with 
a disk. A boundary of length / is specified by a string of / labels q or r; however, these 
labels do not correspond to "spins" in any simple way. Once again, a marked edge on the 
boundary is either attached to a triangle or to another boundary edge. Taking into account 
all possible combinations of labelings for the marked edge and the neighboring triangle or 
boundary edge, we find that the generating function $(g, r) satisfies 



$(g, r) = 1 7(g$r$ + r<l>g$) + a;r<l>r<l> 

+ {a'q + 6'r)(92$ + {b"q + d"r){dgdr^ + 9,9,$) + {d'q + e'r)d^^ 

where we have defined the new variables 

a' = /3a + '-fb , 

b' = ab + 7a b" = (3b + 'yd , 

d' = pd + 76 d" = ad + ■yb , 

e' = ae + 'jd . 



(61) 



(62) 
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The form of these new variables has a geometric interpretation. As an example consider c?', 
the coefficient of qd^^. This term originates from a marked edge labeled q being attached 
to a triangle, which when removed leaves two edges each labeled by r. In such a situation 
the triangle to be removed could have been labeled either qrr or rrr. In the first case the 
weight is given by multiplying the weight d associated with an qrr triangle by the weight /3 
which propagates q to q, while in the second case we multiply e (for a rrr triangle) by 7 (for 
an q-r propagator). It is therefore natural to define d! = (3d + 76, and similarly for the rest 
of (0). 

It is straightforward to check that special cases of the generating equation (^) yield the 
previous equations (0) for the Ising theory and (|i9| ) for the dual theory. 

5.2 Homogeneous boundary conditions 

As with the models studied in the previous sections, there are no general techniques for 
solving the generating equation (|6T|), but we can derive a polynomial equation for the gener- 
ating function restricted to homogeneous boundary conditions (all Q's on the boundary). 
Unfortunately, in the general 2-matrix model the resulting equation for is a quartic with 
extremely complicated coefficients, which it is impractical to present explicitly. Here we will 
describe the steps followed in deriving the quartic. 
We begin by expanding $ as 

$(g, r) = I + Pqq + Prr + Pqqq'^ + Prrr"^ + Pqr{qr + rq) + pgggq^ + Prrrr^ 

(63) 

+Pqqr{qqr + qrq + rqq) + Pqrr{qrr + rqr + rrq) + . . . 

We then define generating functions which depend on q alone, as in the previous sections: 

</'Mg,r)(g) = [w(<9g,9^)$]^^o , (64) 

To solve for 0, we would like to find a set of equations which closes on a set of these 
generating functions. It turns out that there are 10 equations which close on the quantities 
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') fri frri frqri frrri frqqri frqrri frrrr J 



2 + d'(t)rr + a'dlcj) + 2h"dq(t)r 

-iqcj)^ + e'<Prr + h'dl<P + 2d"d, 



'JPq(j) + (3(j)r + (3q(f)dq(f)r + V (prqqr + d' (j). 



b"d, 



qYrr 

+ a'dUr + h"d, 



d, 



qfrr 



■'rqr 



jfqr 



+ b"4>rqrr + d' (, 



a(j) + 7g00r + d"(j)rqr 
Pqa(j) + 70.r + d"(f)rqqr + c' (prqrr + 



rqrr < ^q 
1^1 



qYrqr 



+ a'dlchrr + h"d, 



qfrrr 



+ e'Srrr + h' BlSr + d" 8, 



(65) 



br + h'dl(t)r + d"d, 



qfrqr 



Pra(j) + a(j)r + '~)q(j)(j)rr + d" (prqrr + e'(j)rrrr + b'dg(j)rr + d"dq(j)rrr 



270,+ 

2Q;0r. + 



^ + a'cprqqr + 21)" (j) rqrr + C?'*) 
^ + &'0rggr + 2d" (j)rarr + e'(/ 



Jrqrr 



The first eight of these equations come from considering terms in (|6T|) which begin with a cer- 
tain string; for example the third equation (for dq(j)r) comes from terms of the form qrqqq . . .. 
The last two equations come from terms of the form qrqqq . . . qqqr and rrqqq . . . qqqr, which 
are related by cyclic symmetry to the functions (prqr and (prrr, respectively. 

These equations may be transformed into algebraic relations by replacing the derivatives 
as in the previous two sections; the resulting system of equations can be used to derive a 
single quartic equation involving only 0. 

Note that we can cut down somewhat on the number of undetermined variables by looking 
at equations (p5| ) order by order. These give the following relations between the constants: 

Pq = d'prr + Ct'Pgg + 2h"pqr 

Pqq = P + d'pqrr + Ct'Pggg + 2h" Pqqr 

Pr = e'prr + h'pqq + 2d"pqr 

[bb) 

Pqr = 7 + e'Pqrr + b'Pqgg + 2d"Pqqr 

Pqr T ~l~ 26 Pqrr H" d Prrr H" ^ Pqqr 

Prr ^ ~l~ Pqrr ^ Prrr ~l~ ^ Pqqr 

These allow {pqq,Prr,Pqqq, Pqqr, Pqrr, Prrr) to be expressed in terms of the three constants 
{pq,Pr,Pqr), although we shall not present the explicit results here. 



5.3 Models interpolating between Ising and dual models 

In Sections |] and ^ we discussed models in which Ising spins were coupled directly to 
triangulations and their duals. These two geometrically distinct formulations of the c = 1/2 
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theory are related by a simple coordinate transformation, as was discussed in section ^ 
In this section, we explore a continuous 1-parameter family of models which interpolate 
between the Ising and dual Ising theories. This entire family of theories are related to the 
original Ising theory by linear transformations on the defining hermitian matrices. The 
homogeneous disk amplitude in these interpolating theories can be interpreted physically as 
the disk amplitude in the original Ising theory where an external magnetic field is imposed on 
the boundary. Recently, there has been much work on understanding the effects of boundary 
magnetic fields on the Ising CFT in flat space [^, |61], the discussion here provides a 
formalism by which similar questions can be asked of the theory after coupling to 2D gravity. 
A more detailed investigation of the physical properties of this theory will be carried out in 
a separate publication . 

The original Ising theory is defined by the action 

-1 

S = - Tr(t/^ + V'^) - cTi UV - ^ Tt{U^ + V^) , (67) 

while the dual theory is defined by 

^ = ^^TrX2 + l±^TrF2_^Tr(X3 + 3XF2) . (68) 
2 2 3 



The two models are related by the transformation ([iq ). We now consider the more general 
continuous family of theories defined by the following transformation of the variables of the 
original Ising theory: 

Q = sU + tV , 

(69) 

R = sV + t'V . 

We obtain a theory of the two matrices Q and i?, with action given by ([57|), where the 
parameters are expressed in terms of s,t,s',t' and the couplings c and g of the original 
model: 

a = -^(s" + t" + 2cs't') 
P = -^{s^ + f + 2cst) 
7 = ^(ss' + tt' + est + cs't') 



^ {e - s"') (70) 



/^3' 



(sf - ts 



{sf - ts'] 



{sf - ts'] 



^ 's^ - 1^) 



{sf - ts' 
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The original Ising theory corresponds to s = t' = 1, t = s' = 0, while the dual Ising theory 
corresponds to s = t = —s' = t' = 1/a/2, with the identification g = g/V2. 

The generating equation for this set of models may be derived from the geometric argu- 
ments used in previous sections, or directly from substitution of ( [70D into (|6TD: 

$(g,r) = IH ^\{s^ + + 2cst)q^q^ 

+ {ss' + W + est' + cts')(g$r$ + r$g$) 
+ is''^ + t''^ + 2es't')r^r^ 

ess'^ + ts'^ + st'^ + ctt'2)g9j$ 

+ {s' + t'){es'^ + s't' - es't' + ct'^yd^^ 
-(cs^s' + sts' + stt' + etH')q{dqdr^ + ^^^5$) 
-(css'2 + ss't' + ts't' + ctt'2)r(9g9^$ + 
+ (s + t)(cs2 + - est + et^)qdl^ 
+ {es'^s' + t^s' + sH' + etH')rdl^ 
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{1 - e^){st' - ts'y 



(71) 



In a similar fashion, we may substitute ( ffOD into the quartic satisfied by the generating 
function (p for homogeneous boundary conditions in the general 2-matrix model, to obtain 
a new quartic which is somewhat simpler. Since describes boundary conditions which are 
all Q's, we expect that the quartic will be independent of s' and t' (which merely specify 
R in terms of the original variables U and V^), and this turns out to be the case. We do 
not present the quartic here, but will examine it as part of a discussion of the interpolating 
models in Ref. . 



The generating function 0(g) describes the amplitude for a disk with all Q's on the 
boundary. Since Q = sU + tV, the coefficient of g" in the expansion for is the correlation 
function 

'^TT{sU + tV)A . (72) 



We can therefore calculate the disk amplitude for the original Ising theory with a family of 
different boundary conditions, ranging from homogeneous (s = 1, t = 0) to completely free 
{s = t). Writing 

s = ae^ 

(73) 

t = ae-'' 

where a is a normalization constant, we see that each choice of s,t corresponds to a choice 
of boundary magnetic field. Homogeneous boundary conditions correspond to an infinite 
boundary field, while free boundary conditions correspond to a vanishing boundary field. 
We shall discuss the continuum limit and interpretations of these boundary conditions in 
Ref. 
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6 3-state Potts model 



As mentioned previously, the approach to discretized gravity theories via generating equa- 
tions in free variables, when applied to matrix model theories, may be thought of as an 
extension of the technique of Schwinger-Dyson equations. Although such equations are in 
principle applicable to models where conventional methods (such as orthogonal polynomials) 
are inadequate, in practice it has proven difficult to take advantage of this greater generality. 
One set of models for which orthogonal polynomials seem to fail is those in which the tar- 
get space graph (defining the matter fields) contains a nontrivial cycle, since the nontrivial 
cycle prevents simultaneous diagonalization of the matrices [0, There is some evidence 



that loop equation methods also fail for these models,^ although to our knowledge a clear 
understanding of why this failure occurs is still lacking. 

The simplest theory with a nontrivial cycle is the 3-state Potts model. In fact this model 
is solvable by conventional matrix model techniques, as there exists a change of variables 
which allows the matrices to be diagonalized [^, presumably the same change of 

variables would allow solution by the methods considered here. Nevertheless, in order to 
illustrate the difficulties encountered when a nontrivial cycle is present, in this section we 
examine the 3-state Potts model in its original form, and discuss the obstacles to finding a 
solution for the homogeneous disk amplitude. 

In the 3-state Potts model we associate to each triangle one of three "spins," labelled by 
the free variables Xq, Xi and X2- There is a nearest-neighbor interaction which is the same 
for any two neighboring unequal spins. (The target graph may therefore be thought of as 
a triangle with the three free variables labelling the vertices.) The weight associated with 
each triangulation A is 

^^^^ = , (74) 

where i^(A) is the total number of edges which connect unlike spins. The partition function 
may be cast in matrix model form, in which case it becomes an integral over three matrices 
Xq, Xi, X2, with an action given by 

S = 2(i+c-2c^) Tr [(1 + c)(Xg + Xf + XI) - 2c(XoXi + XqX^ + X.X^)] ^^^^ 
+ ITr{Xi + Xf + XI) . 

Once again we have chosen the quadratic terms in the action such that the propagator takes 
on a simple form. 

From the geometric arguments familiar from previous sections, we derive a generating 



^We would like to thank Matthias Staudacher for explaining this issue to us. 
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equation governing the disk amplitude in this model: 

^{xo,Xi,X2) = 1 + Xo$(a;o + cxi + 0x2)$ + Xi<l>(cxo + Xi + 0x2)$ 

+X2$(cxo + cxi + X2)$ + g{xo + cxi + cx2)<9o$ (76) 
+g{cxo + xi + cx2)dj^ + g{cxo + cxi + X2)9|$ . 

Continuing to follow the procedure outlined in previous sections, from (|76|) we derive a set 
of seven equations in seven functions (0, 0i, 0iii, 0i22, ^iiii, 01122) of the single variable 
Xq: 






= 1 + 2gcxQ(t)ii + gx^dlcj) + xl(t? 


01 


= g{l + c)(f)ii + gcd^(j) + cxo(f)^ 


011 


= 5'c0i22 + 5'0iii + fi'cc?o0i + ca;o00i + 


^001 


= ^-£0122 + fi'c0iii + 5'<9o0i + a;o00i + C0 


0111 


= 5'C01122 + 5'01111 + 5'C<9o011 + 01 + Pl0 + CXO0011 


0122 


= 5'0ii22 + fi'c0iiii + gcd^cpu + C01 + cpi(f) + ca;o00ii 


^0011 


= 5'c0ii22 + 5'c0iiii + gd^cpu + C01 + cpi(f) + a;o00ii 



We might expect to be able to solve these equations to get a single equation for 4>{xo). 
However, closer inspection reveals that they are not independent, with one redundancy 
hidden among them. 

It is possible to uncover a geometric understanding of why these equations are redundant. 
Two of the equations - the one for 0i and the one for 9o0i - can be thought of as governing 
what happens when an edge is removed from a disk with a single Xi and a number of 
Xq's; the first comes from removing the Xi edge, and the second from removing the Xq edge 
immediately next to the Xi. We could imagine acting on such a disk by removing both the 
Xi edge and the xq next to it, in either order; the resulting expressions must be compatible, 
and must be equivalent to performing the appropriate operations on the right hand sides 
of the two equations. But the result of performing such operations on the right hand sides 
is already contained in the other equations in this set (^), so there must be a consistency 
condition, reducing the number of independent equations by one. On the other hand, the 
two equations for (pu and 9o0ii do not share this problem, since the right hand side of 
the equation for 9o0ii contains 0iiii and 0ii22, and the act of removing an edge from the 
triangulations represented by these quantities is not contained in the seven equations. We 
therefore are left with six equations in seven unknowns, insufficient to find a single equation 
for 0. 

By expanding the set of unknowns to include the three additional quantities (0iiiii, 
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</'ii2ii5 <^ii222)5 we can derive four additional equations: 

01111 = 5'C011222 + fi'011111 + 5'C<9o0111 + 011 + Pl01 + Pll0 + 0X000111 
01122 = ^"^011211 + 5'011222 + 5'C<9o0122 + 011 + CPl0l + Cpi20 + CXo00122 

(78) 

<9o0111 = 5'C011222 + fi'C011111 + 5"9o0111 + C0ii + Cpi0i + Cpu(f) + a;o00iii 
^00122 = fi'C011211 + fi'C0ll222 + fi'f^O0122 + C011 + Cpi0i + Cpi20 + a;o00122 

But once again, there is an additional consistency condition which reduces the number of 
independent equations by one. In this case, we have introduced equations with 0iiii and 
01122 on the left hand sides, which leads to another condition on the moves represented by 
the 011 and c?o0ii equations, as discussed above. 

We are therefore still left with one fewer equation than unknowns. As far as we can 
determine, this will continue to be the case no matter how many additional equations we 
introduce, although we have not been able to find a general proof of this fact. It remains 
an open question whether some other method of solution can yield information about this 
model from the generating equation ([76|) . 



7 Multiple Ising spins 

The major limitation on the usefulness of matrix models is the fact that no known methods 
for solving matrix models apply to theories with central charge c > 1. Since the systems 
of the greatest physical interest fall into this category, the "c = 1 wall" has always been a 
major obstacle in the further understanding of 2D gravity and string theories. However, the 
fact that it is difficult to integrate out angular degrees of freedom in matrix models with 
c > 1 has always seemed to be a rather technical problem, and it has been unclear whether 
some other method might make it possible to find analytic results or solutions for systems 
with larger central charge. 

In this section we shall discuss a class of matrix models formed by taking k independent 
Ising spins coupled to quantum gravity. The model with k spins corresponds to a CFT with 
c = k/2 coupled to gravity, and can be described as a 2'^-matrix model. We analyze the 
4-matrix model corresponding to the 2-spin c = 1 CFT coupled to gravity. As was the case 
with the 3-state Potts model, we find that it seems to be impossible to find a closed system 
of equations describing the disk amplitude. 

7.1 Multiple spin model 

The multiple spin Ising theory is defined in a manner precisely analogous to the Ising theory 
in Section |^, except that each triangle carries k independent Ising spins. Thus, in the 
language of triangulations, we sum over all triangulations A where each triangle carries k 
spins (Ti, . . . , CTfc. We set the interaction energy to be identical for each type of spin, so that 
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the weight of a given triangulation A is given by 



W{A) 



s{Ay 



jYX(A)^n(A)^i/(A) 



(79) 



where i^(A) is the total number of adjacent pairs of unequal spins with the same index, and 
c is the interaction between unequal spins (note that we the normalization factor 1 — is 
moved into the matrix model action for convenience in this section). 

The matrix model description of this partition function is given by an integral over 2'^ 
matrices X„ 



.ai-a,,, where CPj G {+1, -1}: 

Z = Y: W{A) = / n DX.,...., exp i-NS) 

A 'J f ^ 



(80) 



{o-l,...,(Tfc} 



with 

S-- 



E 

{o"l v.O-fe} 



2(1 -c 



2\k 



E 



{n,...,Tfc} 



81) 

By the same geometric argument as has been used throughout the paper, we can write a 
generating equation for the disk amplitude as a function of free noncommuting variables 

X, 



which takes the very simple form 



$ = 1 + > c 

CTl,...,(Tfc,ri,...,Tfc 



E 



(82) 



7.2 2-spin model 

For simplicity, we shall restrict attention for the remainder of this section to the 2-spin model 
with c = 1. This 4- matrix model has a continuum limit which can be identified with the 



orbifold c = 1 bosonic string at radius r = 1 [pq , pTf . Labeling the 4 matrices according 

to their binary representations for shorthand {i.e., X3 = X2 = X^^, Xi = X ^, 

Xq = X ), the matrix model action is given by 



^ i,j=0 i=0 



with 



Cij 



1-6 



2\2 



(I -c -C ^\ 

-c 1 — c 

-c 1 — c 

V c2 -c -c 1 y 



^3) 



(84) 
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The generating equation (^) becomes 

+C {{X2 + Xi)^{x3 + Xo)$ + (X3 + Xo)^{x2 + Xi)^) 
{x^^Xq^ + X2^Xi^ + Xi^X2^ + Xo^Xs^) 

(85) 

+g {xsdi^ + X2di^ + xidf^ + xod^^) 

+gc ((X3 + xo)idl<i> + 92$) + {x2 + xi){dl<l> + 52$)) 

+gc^ {x3d^^ + X2dl^ + xi(9^$ + xo9|$) . 

As for the 3-state Potts model, it appears not to be possible to derive a single algebraic 
equation for the homogeneous disk amplitude from this equation. Once again, however, we 
do not have any proof that this must be the case. However, it seems unlikely that any 
algebraic equation for the homogeneous disk amplitude exists in view of the logarithmic 
behavior of other previously studied c = 1 models. 

It is nevertheless conceivable that some other means of analyzing the generating equation 
(^) could lead to a single equation characterizing 0. Furthermore, we mention the promising 
result that all coefficients Ps can be determined in terms of the three coefficients po, pooo, 
Poooo; just as in the Ising theory where only two such coefficients were necessary. Presumably, 
however, any method for solving the c = 1 model would not generalize to the c = k/2 models 
with k > 2 since the number of unknowns in the higher k models grows much faster than 
when k = 2. 



8 Dually weighted models 

So far in this paper, all of the specific models which we have considered have simple descrip- 
tions in terms of the matrix model formalism (although we described a class of boundary 
conditions for the Ising theory in Section |3.4| which cannot easily be described using matrix 



model methods) . However, the methodology we have developed is applicable to a much 
wider class of discrete gravity models. In general, we can construct discrete matter fields 
by endowing a dynamical triangulation with an arbitrary set of local discrete data. For 
example, we can simultaneously give discrete labels (/j, Cj, Vi) to faces (polygons), edges, 
and vertices of a triangulation. If we make the weight function local, in the sense that the 
only interactions are between adjacent faces, edges, and vertices, then the partition function 
corresponding to such a theory is of the form 



n exp{p{f,{eedf},{vedfm 

/Gfaces 

n exp{r]{{f ■.eedf},e,{vedem (86) 

eg edges 
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n exp{u{{f -.v edf},{e:v ede},v))) 



fSvertices 



where A is a set of combinatorial data {fi,ei,Vi) on the simphces in a triangulation A; p 
is a function of the geometry and labels of a fixed face / and the edges and vertices on its 
boundary; similarly rj is a function of a fixed edge e, the faces of which e is a boundary, and 
the vertices on the boundary of v; and u is a. function of a vertex v and the edges and faces 
of which it is a boundary. 

A generating equation can be derived for a model of the general form (^) by associating 
independent noncommuting variables with the discrete labels on edges, vertices, and faces. 
Although the complexity of such an equation will increase significantly as the set of discrete 
matter fields increases, in principle the methods of this paper can provide an approach to 
the analysis of any such theory. 

A particular set of theories in this more general class are the "dually weighted graph" 



models recently studied by Di Francesco and Itzykson and by Kazakov, Staudacher, and 



Wynter |3^. In these models, weights are associated with vertices of a triangulation accord- 



ing to the coordination number, and weights are associated with the faces of a triangulation 
according to the number of incident edges as in a general one-matrix model. The dually 
weighted graph models are of particular interest because they contain continuous families of 
theories which interpolate between conformal field theories on flat space and conformal field 
theories completely coupled to Liouville gravity. Although it is possible to construct most 
dually weighted graph models as matrix models by introducing external matrices encoding 
the vertex weights, the resulting models are rather unwieldy and difficult to solve by stan- 
dard matrix model techniques. An example of a model which can be analyzed using matrix 
model methods is the model where vertices of both the original and dual graph are restricted 
to have even coordination number |^ . 



We now derive the generating equation for a general dually weighted graph model. The 
most general dually weighted graph model has a weight of gi for every face with i edges, and 
a weight of Wi for every vertex with coordination number i. Thus, the partition function for 
a general dually weighted graph model is given by 



i7) 



where Fi{A) is the number of z-gon faces in A and Vi(A) is the number of vertices with 
coordination number i. 

To describe the generating function for the disk amplitude of this model, we must as- 
sociate noncommuting variables with the vertices on the boundary of the disk. These non- 
commuting variables must contain information about the coordination number of the vertex, 
as well as the number of faces inside the disk which are incident on that vertex. Thus, we 
choose variables Xi\j, with i < j to denote a vertex with i interior faces and a coordination 
number of j. (An example is shown in Figure |^.) 
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Figure 8: Variables Xi\j on the vertices around the boundary of a disk 
in a dually weighted graph model. The subscript j indicates the coor- 
dination number of the vertex; i indicates the number of faces incident 
on the vertex on the interior of the disk. 



The generating equation for the general duaUy weighted graph model is 

$ = 1 + gkXi+i\jXi>+iij'di'\fd^~^diij^ 

\i<j,i'<j',k 

+ m Xi\kXi>\k'{di>+j>-k'\k'^)xj'ik'Xj\k{di+j-k\k^) (89) 

\i+j>k,i'+j'>k' J 

where we have defined 

i>2 

As an extremely simple example, let us consider a "flat space" dually weighted graph 
model where only hexagonal faces and vertices with coordination number 3 are allowed 
{i.e., gi = 6i^G; Wi = Si^s). This is simply a theory of "triangulations" of flat surfaces by 
hexagons; each term in the disk amplitude for this theory corresponds to a closed path on 
the honeycomb lattice which bounds an embedded oriented surface composed of hexagons. 
The generating equation for this model reduces to (writing Xi — x^^) 

$ = 1 + a;2X2(?i(?2(?i<l> + a;2Xi9o929i$ + a;iX29i929o$ + XiXi^o^g^o^ 



+X2 [a;2(<9i$)a;2 + xi{do^)x2 + X2(<9o$)xi] X2(9i$) 

+X2 [X2{di^)x2 + Xi{do^)x2 + X2{do^)xi] Xi(9o$) 
+Xi [X2{di^)x2 + Xi{do^)x2 + X2ido^)Xi] X2(<9o$). 
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There are a number of slightly more interesting models which are very simple to describe, 
and which connect flat space theories such as the above with fluctuating 2D gravity theories 
which presumably correspond to coupling in some form of matter fields. For example, con- 
sider the theory with only square faces, and vertices of coordination number between 3 and 
5 defined by 

9i = ^i,4 

Wi = (5i,4 + a5i,3 + b6i^5- (91) 

When a = b = 0, this is simply another fiat space theory. As a and b become nonzero, fluc- 
tuations are introduced. It should be possible by tuning the ratio a/6 to find a 1-parameter 
family of theories which correspond to a continuous limit; any detailed understanding of the 
behavior of this class of theories would be extremely interesting for understanding the effects 
of coupling dynamical gravity to a flat theory. The generating equation for this model is 
simply the restriction of (^) to the theory with weights (pi]). 

It would be interesting to see whether for these simple models there is any reasonable 
algebraic method for extracting information from the generating equation. 



9 Conclusion 

We have developed a systematic method for analyzing discretized theories of 2D gravity 
with matter. Our method combines a geometric/combinatorial approach with the tools of 
free variables and loop equations which have previously been used to study matrix mod- 
els. Although this approach gives rise to equations which are essentially equivalent to the 
Schwinger-Dyson loop equations for standard matrix model theories, our results are more 
general in several ways. First, because our methodology is completely based on a geometric 
definition of the discrete gravity theories, the analysis applies to models which cannot be 
described in terms of the matrix model language. Second, we can describe correlation func- 
tions in standard matrix model theories using variables which are apparently inaccessible 



using the matrix model formalism. For example, in |^ we describe the calculation using 
our techniques of the exact correlation function between disorder operators in the Ising the- 
ory. As far as we have been able to ascertain, this correlation function cannot be computed 
in any simple way using matrix model methods. Finally, we believe that the principal result 
presented in this paper, the general single generating equation for the generating function 
of all disk amplitudes, represents a conceptual and algebraic simplification and synthesis of 
previous work in this direction. 

The methods we have presented here have enabled us to analyze in a systematic fashion a 
variety of models describing 2D gravity with matter fields. One question which has never been 
fully resolved is the question of which matrix models allow solution by loop equation methods. 
To date, the only theories which have proven tractable using this approach have been the pure 
gravity and c = 1/2 models. We have found that it is possible to use this approach to derive 
polynomial equations satisfied by the generating function of disk amplitudes not only in the 
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standard Ising formulation of the c = 1/2 theory, but also in the dual formulation, and even 
in a continuous family of models connecting these two formulations. In addition, we have 
found that a very general class of two-matrix models, which cannot be solved by orthogonal 
polynomials, have uniform boundary condition disk amplitude generating functions which 
satisfy a quartic. 

We have attempted to solve several other models using these methods, including the 
3-state Potts model and the c = 1 model defined by two Ising spins coupled to gravity. We 
have found that these models resist solution, in accordance with the empirical belief that 
models whose target graphs contain loops somehow have more intrinsic degrees of freedom, 
and are not soluble by loop equations or algebraic means. However, in our attempt to 
solve these equations, we uncovered algebraic structure to the obstruction, characterized by 
degeneracies in subsets of the loop equations. This algebraic structure may give some hint 
of how to use other approaches to solve these models. 

In the discrete theories that we have considered, the existence of a single equation describ- 
ing the effect of all possible moves that reduce a triangulation suggests that this equation 
should be related to a Hamiltonian operator for the theory. In this context it would be inter- 
esting to consider the continuum limit of the equation in terms of noncommuting variables 
and see if can be related to recent attempts to construct string field hamiltonians for theories 
with c < 1 [22-30]. 

The approach we have detailed here simplifies the calculation of correlation functions 
in a variety of interesting 2D theories. Some particular results along these lines will be 
described in [Q, However, the major conceptual issue which needs to be dealt with for 
these methods to be applicable in a truly general context is the development of some more 
sophisticated and powerful methods for dealing with algebraic equations in noncommuting 
variables such as (|30D . This is a very difficult problem, as has been emphasized in a number 
of recent papers ^ ^ ; however, in these works and others some progress has been 

made, giving hope that in the future these equations may seem more tractable. 
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